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Supercavities  in  Compressible  Fluid 


A.  D.  Vasin 

State  SRC  TsAGI,  Moscow,  Russia 
17,  Radio  str.,  Moscow,  Russia 


Summary 

The  basic  results  obtained  by  the  author  in  the  theory  of  supercavities  in  compressible  fluid  have  been 
presented.  Some  peculiarities  of  the  theoretical  methods  application  for  the  supersonic  flows  and  their 
difference  from  the  analogous  flows  in  air  have  been  noted.  The  basic  correlations  for  the  shock  waves 
in  the  water  are  considered. 

Axissymetric  supercavities  in  sub-  and  supersonic  compressible  fluid  have  been  researched  using  the 
asymptotic  method.  An  integro-differential  equation  for  the  cavity  profile  was  deduced.  The  outer 
solution  was  found  and  for  the  thin  cone-cavitators  the  solution  was  obtained  for  the  whole  area  of  the 
flow. 

Based  of  the  modem  numerical  methods  the  algorithm  of  calculation  of  axisymmetric  cavities  in  sub-  and 
supersonic  compressible  fluid  flow  has  been  worked  out.  The  results  of  calculation  of  axysimmetric 
supercavities  downstream  of  a disk  within  the  range  of  the  Mach  numbers  0<M<1.4  have  been 
represented.  As  a result  of  the  analysis  for  the  different  Mach  numbers  the  following  was  defined:  the 
cavitation  number,  the  cavitation  drag  ratio,  the  mid-section  radius,  the  cavity  shape,  the  distance  of 
shock  wave  from  a disk  (for  supersonic  flow).  The  conical  flows  in  a supersonic  water  stream  have  been 
calculated  by  the  numerical  method. 

The  results  obtained  by  the  author  have  been  compared  with  the  known  experimental  and  theoretical 
results  of  other  investigations. 


Introduction 

The  axisymmetric  and  three-dimensional  supercavities  in  incompressible  fluid  have  been 
intensively  studied  since  the  early  1950’s.  A large  amount  of  works  have  been  published  on  three- 
dimensional  supercavities  in  incompressible  fluid.  There  were  about  several  hundred  of  them.  It  should 
be  noted  that  there  have  been  a few  publications  on  the  three-dimensional  supercavities  in  a compressible 
fluid  (in  contrast  to  the  case  of  incompressible  fluid).  This  situation  is  a result  of  both  the  absence  of 
practical  interest  to  this  problem  before  1980  and  the  complexity  of  the  problem  itself.  The  only  work 
was  Gurevich’s  work  [1]  that  was  published  in  1947.  In  this  work  the  expansion  law  for  a cavity 
downstream  of  a body  was  examined  for  zero  cavitation  number  and  generalized  to  include  the  subsonic 
case. 


Since  the  early  1980’s  in  the  USSR  and  in  the  USA  the  interest  in  supercavities  in  a compressible 
fluid  has  increased.  The  necessity  to  come  up  with  the  supercavities  theory  in  a compressible  fluid 
appeared.  The  following  works  have  been  made  based  on  the  slender  body  theory.  These  are  the  works: 
Yakimov  [2,3],  the  author  [4,5]  and  Serebryakov  [6,7].  On  the  basis  of  the  early  numerical  method  Al’ev 
[8]  has  calculated  the  cavity  past  a thin  cone.  The  theoretical  investigations  on  two-dimensional 
supercavities  in  compressible  fluid  were  carried  out  by  Nishiyama  [9]. 

The  slender  body  theory  has  a limited  application  for  calculation  of  supercavities  in  compressible 
fluid.  For  blunt  cavitators  such  as  a dusk  or  a blunt  cone,  this  theory  does  not  provide  a way  to  obtain  the 
complete  data  for  the  cavitation  flow,  for  example,  the  cavity  shape  immediately  downstream  of  the 
cavitator,  and  the  cavitation  drag.  The  linear  equation  for  the  potential  of  compressible  fluid  flow  past 


Paper  presented  at  the  RTOAVT Lecture  Series  on  “Supercavitating  Flows  ”,  held  at  the  von  Karman 
Institute  (VKI)  in  Brussels,  Belgium,  12-16  February  2001,  and  published  in  RTO  EN-010. 


16-2 


slender  axisymmetric  bodies  is  inapplicable  on  the  transonic  velocity  range.  Moreover,  for  the  case  of 
supersonic  flow  the  slender  body  theory  faces  some  major  difficulties  which  is  shown  by  the  author’s 
investigations. 

The  recent  development  of  the  numerical  methods  and  widespread  use  of  computers  have  made  it 
possible  to  calculate  compressible  cavitation  flows  efficiently.  It  should  be  noted  that  the  numerical 
calculations  have  been  carried  out  by  Zigangareeva  and  Kiselev  [10]  and  by  the  author  [11,12]. 


1.  The  application  of  the  slender  body  theory  for  the  investigation  of  supercavities 
in  subsonic  flow 

The  case  of  subsonic  cavitation  flow  past  a thin  cone  has  been  considered  in  chapter  2 of  lecture 
#1.  The  Riabouchinsky  scheme  has  been  applied  (Fig.  3 in  lecture  #1).  An  integro-differential  equation 
for  the  cavity  profile  has  the  form  of  correlation  (2.6)  in  lecture  #1.  In  lecture  #1  we  have  considered  the 
case  when  the  parameter  of  the  cavity  thinness  has  the  same  order  as  the  parameter  of  the  cavitator 
thinness  £ equal  to  the  ratio  of  the  cone  radius  to  its  height.  For  the  whole  area  of  the  flow  the  solution 
has  been  obtained  by  the  method  of  expansion  into  asymptotic  rows  (rows  have  the  form  of  correlations 
(2.6)  in  lecture  #1).  The  first  approximation  of  the  cavity  profile  is  represented  by  an  ellipsoid  of 
revolution  (equation  (2.9)  in  lecture  #1).  The  second  approximation  is  determined  from  the  equation 
which  was  obtained  after  the  equalization  of  members  at  e4  [13]. 


Fig.  1 Cavity  profile  in  subsonic  flow  for  the  cone  oc= 5°;  the  broken  curve  corresponds  to  the 
first  approximation,  the  continuous  curve  corresponds  to  the  second  approximation 


In  Fig.  I the  broken  curve  corresponds  to  the  first  approximation  of  the  cavity  profile  R(x)  past 
the  cone  with  the  apex  half-angle  equal  to  5°.  The  total  length  L equals  200  (the  geometric  dimensions 
are  scaled  by  the  radius  of  the  cone  base,  the  description  of  the  geometric  dimensions  is  given  in  chapter 
2 of  lecture  #1).  The  continuous  curve  corresponds  to  the  second  approximation  of  the  cavity  profile. 
From  Fig.  1 it  is  evident  that  the  second  approximation  is  close  to  the  first  one,  i.e.  the  second  members  of 
asymptotic  row  are  smaller  than  the  first  ones. 

Besides,  for  the  cavities  past  the  thin  cones  the  cavitation  numbers  <J  are  determined  at  the 
different  Mach  numbers  M (<7  = 2 (P0  - Pk)/ pfjj2,  where  Pu  is  the  free  stream  pressure.  PK  is  the 

pressure  within  the  cavity,  p is  the  fluid  density,  V0  is  the  free  stream  velocity;  M=V(/ao,  where  a0  is  the 
free  stream  speed  of  sound).  The  dependences  o(M)  are  determined  at  the  constant  cavity  length  Lk  and 
the  total  length  L equal  to  200. 
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In  Fig.2  the  dependences  o(M)  are  shown  at  various  cone  apex  half-angles  a (curves  1-3 
correspond  to  a=15°,  10°  and  5°  respectively).  Thus,  for  the  given  cone  at  the  constant  value  of  L or  Lk 
the  cavitation  number  increases  with  the  Mach  number  and  the  cavity  profile  remains  invariable. 

After  implementation  of  the  transformations  analogous  to  the  ones  in  chapter  2 of  lecture  1 we 
can  determine  the  expression  for  the  pressure  coefficient  on  the  thin  cone  surface  [13].  Integrating  the 
pressure  coefficient  on  the  cone  surface  we  obtain  the  cavitation  drag  coefficient  Cx  (normalized  by  the 
area  of  the  cone  base  kR]  , where  R„  is  the  radius  of  the  cone  base).  In  Fig.  3 the  dependences  CX(M)  are 
shown  at  various  cone  apex  half-angles  (curves  1-3  correspond  to  a=15°,  10°  and  5°  respectively). 


Fig.  3 Dependence  Cx  on  M;  curves  1-3  correspond  to  a=15°,  10°,  5' 
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The  character  of  the  dependences  CX(M)  (Fig.3)  corresponds  to  one  of  the  dependences  a(M) 
(Fig.2).  This  is  a consequence  of  the  law  of  conservation  of  momentum  applied  to  the  cavity  middle- 
section  plane.  The  law  of  conservation  of  momentum  can  be  written  in  the  following  form  [14] 


R 


k 


(1.1) 


where  RK  is  the  radius  of  the  cavity  mid-section  normalized  by  the  radius  of  the  cone  base,  k is 
the  coefficient  which  value  is  close  to  1 at  the  small  cavitation  numbers.  In  the  first  approximation  we 
can  neglect  the  dependence  of  k on  the  Mach  number.  The  solution  for  the  thin  cones  has  shown  that  at 
the  constant  value  of  Lk  or  L the  value  of  RK  remains  the  same  at  the  variation  of  the  Mach  number.  Then 
according  to  formula  (1.1)  the  dependences  CX(M)  and  a(M)  must  satisfy  the  following  correlation 


cJm) 

o(M ) 


= const 


(1.2) 


where  the  constant  is  the  function  of  a and  Lk.  The  calculations  have  shown  that  correlation  (1.2)  is 
satisfied  with  sufficient  accuracy  at  the  Mach  numbers  on  the  interval  0<M<0.8  [13],  The  cavity  past  a 
thin  cone  in  incompressible  fluid  (M=  0)  has  been  considered  as  a particular  case.  The  value  of  the  cone 
apex  half-angle  equals  15°.  For  this  cone  the  results  obtained  based  on  the  slender  body  theory  have  been 
compared  with  the  results  of  the  numerical  calculation  [15]  for  incompressible  fluid.  The  comparison  has 
shown  good  agreement  of  these  results  [13], 

Let  us  consider  the  case  when  the  parameter  of  the  cavity  thinness  is  far  less  than  the  parameter 
of  the  cavitator  thinness.  In  this  case  we  neglect  the  cavitator  dimensions  in  comparison  with  the  cavity 
dimensions.  In  order  to  determine  the  outer  solution  we  assume  that  the  radius  of  the  cavitator  equals 
zero.  Let  us  place  the  origin  of  coordinates  in  the  middle  of  the  cavity  length  and  the  geometric 
dimensions  normalize  by  the  semi-length  Lk/2.  In  this  case  integro-differential  equation  (2.6)  in  chapter  2 
of  lecture  #1  takes  the  fonn  [4,13] 


ti 


2 u 


■ + u In 


„ (l-M2)w  r u"{xx  )-  u"(x)  l u'(- 1)  u'( l) 


4(1  — *!) 


-1 


x- 


■dx]  - — — — + — = 2a 
1 + x 1 - x 


R2,  R(±  l)=0,  R{ 0)=£, 


(1.3) 


where  £,  = 2 Rk  / Lk  - 1 / A is  the  small  parameter  equal  to  the  inverse  amount  of  the  cavity 
aspect  ratio  X. 

We  seek  the  outer  solution  in  the  fonn  of  expansion  into  asymptotic  rows  with  the  small 
parameter  E\.  For  the  squared  cavity  radius  and  the  cavitation  number  these  rows  have  the  form 
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After  substitution  of  (1.4)  in  (1.3)  and  conservation  of  the  two  members  of  the  rows  integro- 
differential  equation  (1.3)  is  transformed  to  two  differential  equations.  Solving  the  differential  equations 
we  obtain  the  expression  for  the  cavity  profile  and  the  dependence  of  the  cavitation  number  on  the  cavity 
aspect  ratio  and  the  Mach  number  [4,13] 
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where  e is  the  base  of  natural  logarithms. 


(1.5) 


The  first  approximation  of  the  cavity  profile  is  an  ellipsoid  of  revolution  (as  in  incompressible 
fluid).  From  the  first  expression  (1.5)  it  follows  that  the  added  member  to  an  ellipsoid  is  small  at  the 
small  cavitation  numbers,  for  example,  then  A=14.66  and  cr=0. 02  (M= 0)  the  difference  between  the 
cavity  shape  and  an  ellipsoid  of  revolution  does  not  exceed  1%.  In  general  the  cavity  aspect  ratio  is 
dependent  upon  the  cavitation  and  Mach  numbers.  From  the  second  expression  (1.5)  it  follows  that  at  the 
constant  value  of  the  cavity  aspect  ratio  the  cavitation  number  increases  along  with  the  Mach  number. 
The  analogous  dependence  is  obtained  in  the  case  of  the  subsonic  cavitation  flow  past  a thin  cone.  As  the 
cavitation  number  decreases  the  influence  of  fluid  compressibility  on  the  cavity  aspect  ratio  becomes 
insignificant. 

Thus,  the  investigation  carried  out  on  the  basis  of  the  slender  body  theory  has  shown  that  in 
subsonic  flow  we  can  neglect  the  influence  of  the  Mach  number  on  the  cavity  shape:  the  cavity  shape  is 
described  by  an  ellipsoid  of  revolution  with  sufficient  accuracy  (as  in  incompressible  fluid).  An  ellipsoid 
of  revolution  is  obtained  from  the  outer  solution  (1.5)  (when  we  neglect  the  cavitator  dimensions  in 
comparison  with  the  cavity  dimensions)  and  from  the  solution  for  the  thin  cones  ((2.9)  in  lecture  #1) 
(when  the  parameter  of  the  cavity  thinness  has  the  same  order  as  the  parameter  of  the  cavitator  thinness). 
The  Mach  number  mainly  influences  the  value  of  the  cavitation  drag  of  the  cavitator  (the  cavitation  drag 
coefficient  increases  along  with  the  Mach  number).  At  the  constant  cavity  length  the  cavitation  number 
increases  along  with  the  Mach  number  according  to  the  law  of  conservation  of  momentum. 


2.  The  application  of  the  slender  body  theory  for  the  investigation  of  supercavities 
in  supersonic  flow 

Let  us  consider  the  cavitation  flow  past  a thin  cone  in  supersonic  stream.  We  apply  the 
Riabouchinsky  scheme  (Fig.  3 in  lecture  #1).  The  linearized  equation  for  the  perturbed  velocity  potential 
in  the  cylindrical  coordinate  system  is  written  as  following  [16] 
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where  (p  = <jO  * / V0Rn  is  the  dimensionless  perturbed  velocity  potential,  (p*  is  the  perturbed 
velocity  potential. 

The  geometric  dimensions  are  scaled  by  the  radius  of  the  cone  base  (the  description  of  the 
geometric  dimensions  is  given  in  chapter  2 of  lecture  #1).  For  supersonic  flow  the  perturbed  velocity 
potential  that  satisfies  equation  (2.1)  can  be  written  in  the  form  of  integral  [16] 


<P  = 


. 1-lS^r  / \ , 

J_  r 
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(2.2) 


where  q(x \)  is  the  intensity  of  sources  and  sinks  on  the  axis  of  symmetry.  Near  the  surfaces  of 
the  thin  cavitator  and  the  cavity  potential  (2.2)  that  satisfies  the  kinematic  boundary  condition  has  the 
following  asymptotic  representation  [ 1 7] 


4/r  4x 
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(2.3) 


where  S'(x ) is  the  derivative  of  the  dimensionless  area  of  the  cross  section  with  respect  to  x 
coordinate. 

For  supersonic  flow  the  dynamical  boundary  condition  on  the  cavity  surface  has  the  same  form  as 
for  subsonic  flow  [5]  and  is  written  as  following 


9x  [ dr 
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Substituting  (2.3)  in  (2.4)  we  obtain  the  integro-differential  equation  for  the  cavity  profile  past  a 
thin  cone  in  supersonic  flow 
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We  assume  that  the  cone  has  the  same  thinness  order  as  the  cavity.  The  thinness  order  of  the 
cone  e equals  the  ratio  of  the  radius  of  the  cone  base  to  its  height.  We  seek  the  solution  for  the  whole 
area  of  the  flow  in  the  form  of  expansion  into  asymptotic  rows  with  the  small  parameter  e.  The  rows 
have  the  form  of  correlations  (1.4)  if  we  substitute  e with  £\.  The  boundary  conditions  are  written  in  the 
same  form  as  for  subsonic  flow  [7,13] 

x-l : R = 1,  R'  = e 

x = l + Lk  :R  = 1 
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After  substitution  of  the  asymptotic  rows  in  (2.5)  and  conservation  of  the  two  members  of  rows 
integro-differential  equation  (2.5)  is  transformed  to  two  differential  equations.  The  first  differential 
equation  with  the  boundary  conditions  is  written  as  follows 


d2R20 
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(2.6) 


The  solution  of  equation  (2.6)  has  the  form 


R 2 = <7,  (a  - x)(x  -b) 


From  expression  (2.7)  it  follows  that  an  ellipsoid  of  revolution  is  the  first  approximation  of  the 
cavity  shape.  The  analogous  result  has  been  obtained  in  subsonic  flow  (expression  (2.9)  in  lecture  #1). 
The  second  approximation  is  determined  from  the  differential  equation  which  was  obtained  after 
equalization  of  members  at  e4. 

In  Fig.  4 the  broken  curve  represents  the  first  approximation  of  the  cavity  profile  past  the  cone 
with  the  apex  half-angle  equal  to  5°.  The  total  length  equals  200.  The  continuous  curve  represents  the 
second  approximation  of  the  cavity  profile. 


8. 
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Fig.  4 Cavity  profile  in  supersonic  flow  for  the  cone  a=5°;  the  broken  curve  corresponds  to  the 
first  approximation,  the  continuous  curve  corresponds  to  the  second  approximation 


Figure  4 is  analogous  to  Fig.  1 that  illustrates  the  cavity  profile  past  the  same  cone  in  subsonic 
flow.  From  Fig.  4 it  follows  that  as  distinct  from  subsonic  flow  (Fig.  1),  in  supersonic  flow  the  second 
approximation  of  the  cavity  profile  essentially  differs  from  the  first  one,  and  the  second  members  of 
asymptotic  row  are  not  smaller  than  the  first  ones.  Besides,  the  unreal  shape  of  the  cavity  contradicts 
with  the  physical  laws  of  cavitation  flows. 

The  analysis  of  this  effect  has  shown  that  the  slender  body  theory  is  not  valid  for  the  flow  past 
the  edge  of  the  cone.  In  the  supersonic  flow  the  Prandtl-Meyer  stream  arises  at  the  edge.  The  stream 
must  turn  at  the  angle  that  is  dependent  on  the  pressure  within  the  cavity  or  the  cavitation  number.  The 


16-8 


incline  of  a free  stream  line  is  less  than  the  cone  apex  half-angle.  It  is  impossible  to  satisfy  this  condition 
within  the  range  of  the  slender  body  theory.  The  first  peculiarity  of  the  slender  body  theory  application 
to  the  supersonic  cavitation  flow  is  its  discrepancy  with  the  boundary  condition  at  the  cavitator  edge. 

In  order  to  determine  the  solution  it  is  necessary  to  consider  the  supersonic  flow  around  the  cone, 
obtain  the  Mach  number  at  its  edge  and  match  the  Prandtl-Meyer  stream  with  the  cavitation  flow.  It  is 
very  difficult  to  carry  out,  as  the  cavity  shape  and  the  cavitation  number  are  not  known  beforehand. 

Let  us  determine  the  pressure  coefficient  on  the  surface  of  the  thin  cone  in  supersonic  flow.  The 
expression  for  the  pressure  coefficient  has  the  form  [ 1 6] 


p dx  dr  , 


(2.8) 


where  C = l{Pcon  - P0 )/ pV02  is  the  pressure  coefficient,  Pco„  is  the  pressure  on  the 


cone 


surface. 


After  differentiating  expression  (2.3)  the  components  of  the  perturbed  velocity  on  the  cone 
surface  are  written  as  following 


d(p  e2  |n(M2  -\)e 
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(2.9) 


After  substitution  of  (2.9)  in  (2.8)  for  the  pressure  coefficient  we  obtain  the  following  expression 


ln(M2  -l)g2 
4 


(2.10) 


Formula  (2.10)  coincides  with  the  expression  that  defines  the  pressure  coefficient  on  the  thin 
cone  in  supersonic  flow  of  air  [16].  The  pressure  coefficient  is  connected  to  the  cavitation  drag 
coefficient  by  the  following  dependence 


i 

Cx  = 2eJ  CpR\  (x)clx  + a 


(2.11) 


where  Cx  is  the  cavitation  drag  coefficient  normalized  by  the  area  of  the  cone  base;  /(,  (x)=  £X 
is  the  radius  of  the  cone  cross  section. 

The  pressure  coefficient  is  constant  on  the  cone  surface  and  expression  (2.11)  takes  the  form 


Cx=Cp+a 


(2.12) 
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We  can  consider  the  coefficient  Cp  to  be  the  cavitation  drag  coefficient  at  the  zero  cavitation 
number  Cxo.  The  formula  is  written  as  follows 


C,  = C + o 


(2.13) 


Expression  (2.13)  is  the  dependence  of  the  cavitation  drag  coefficient  on  the  cavitation  number  at 
the  constant  Mach  number. 

Let  us  consider  the  case  when  the  parameter  of  the  cavity  thinness  is  far  less  than  the  parameter 
of  the  cavitator  thinness.  A similar  case  was  considered  for  subsonic  flow  in  chapter  1.  Let  us  place  the 
origin  of  coordinates  in  the  middle  of  the  cavity  length  and  normalize  the  geometric  dimensions  by  the 
semi-length  L/l.  Integra-differential  equation  (2.5)  takes  the  form 


u'2  u\  (m2  -l)« 

+ — In- — 

4 u 2 4(x  + l) 


r u"(x[  ) - u"(x) 


-dx. 


-1 


u = R2,  R{±  l)=0,  R( 0)  = £, 


x + 1 


-o 


(2.14) 


where  £,  = 2 Rk  / Lk=\  i A is  the  small  parameter  equal  to  the  inverse  amount  of  the  cavity 
aspect  ratio  A. 

For  determination  of  the  outer  solution  we  neglect  the  cavitator  dimensions  in  comparison  with 
the  cavity  dimensions  (the  radius  of  the  caivtator  equals  zero  in  (2.14)).  We  seek  the  outer  solution  in  the 
form  of  expansion  into  asymptotic  rows  with  the  small  parameter  E\.  For  the  cavity  radius  squared  and 
the  cavitation  number  these  rows  have  the  form  (1.4).  After  substitution  of  (1.4)  in  (2.14)  and 
conservation  of  the  two  members  of  the  rows  integro-differential  equation  (2.14)  is  transformed  to  two 
differential  equations.  Solving  the  differential  equations  we  obtain  the  expression  for  the  cavity  profile 
and  the  cavitation  number  dependence  on  the  cavity  aspect  ratio  and  the  Mach  number 


R2  — £,2 


(i-*4 


+ - 


x2  In  4 + In 


(l  + xf-^^l-x)^2)]] 


2 In  A 


<7  = — — In 


A 


Vca/m2  - 1 


(2.15) 


The  comparison  of  the  dependences  (2.15)  with  the  analogous  ones  (1.5)  obtained  for  subsonic 
flow  shows  that  in  supersonic  flow  the  first  approximation  of  the  cavity  shape  is  an  ellipsoid  of 
revolution  (as  in  incompressible  fluid).  The  numerator  of  the  second  member  in  the  asymptotic  row 
(2.15)  can  be  presented  as  the  sum  of  two  functions  f(x)  and  f2(x) 


f\ (x)=  x2  In  4 - ln[(l  + x)^1+'^(l  - x)(l  1 J 
f2  (x)=  In  (l  + xf  _1^(l  - x)*'!  ’ ' 


The  function//(x)  corresponds  to  the  second  approximation  of  the  cavity  profile  in  subsonic  flow 
(1.5).  In  case  of  supersonic  flow  the  odd  function  f2(x)  is  added  to  the  even  function  f(x)  and  the  cavity 
becomes  asymmetric  relatively  to  the  middle  of  its  length.  The  cavity  asymmetry  is  explained  by  the  fact 
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that  in  supersonic  flow  the  influence  of  sources  and  sinks  spreads  downstream  only.  The  expression  for 
f2(x)  shows  that/2(xj>0  at  x<0  and  f2(x)< 0 atx>0,  i.e.  the  cavity  is  thicker  in  the  head  than  in  the  back  part 
and  the  mid-section  displaces  ahead. 

The  calculations  have  shown  that  at  the  cavitation  numbers  that  have  the  order  10'2  the  second 
approximation  of  the  cavity  shape  in  (2.15)  is  close  to  the  first  one  and  the  difference  does  not  exceed 
several  percent.  From  the  second  expression  (2. 1 5)  it  follows  that  the  cavity  aspect  ratio  is  dependent 
upon  the  cavitation  and  Mach  numbers.  As  the  cavitation  number  decreases  the  influence  of  the  fluid 
compressibility  on  the  cavity  aspect  ratio  becomes  insignificant. 


3.  The  application  of  the  numerical  methods  for  calculation  of  supercavities  in 
subsonic  flow 

The  modem  numerical  methods  enable  us  to  calculate  the  cavitation  flows  in  compressible  fluid 
efficiently.  The  author  has  worked  out  the  algorithm  of  calculation  of  axisymmetric  cavities  in  sub-  and 
supersonic  compressible  fluid  flow  over  a wide  range  of  cavitation  numbers  [11,12], 

Let  us  consider  axisymmetric  subsonic  cavitation  flow  past  a disk  in  accordance  with  the 
Riabouchinsky  scheme;  in  this  scheme  the  cavity  is  closed  by  a disk  of  the  same  dimensions  as  the 
cavitator  disk  (Fig.  5) 


Fig.  5 Flow  schematic  and  computation  domain.  1 is  the  forward  shock,  2 is  the  trailing  shock 


In  view  of  the  flow  symmetry,  we  will  examine  the  flow  in  a meridional  plane.  The  problem 
consist  in  determining  the  shape  of  the  cavity  boundary  MN  that  satisfy  the  constant  velocity  and 
impermeability  conditions  for  a given  cavity  length  Lk  (see  Fig.  5 where  AM  is  the  disk-cavitator,  NB  is 
the  closing  disk,  CFED  is  the  external  boundary).  The  geometric  dimensions  are  scaled  by  the  disk 
radius.  The  external  boundary  CFED  consists  of  three  parts:  CF  and  ED  are  parts  of  the  circles  with 
radius  120  and  centers  at  the  points  A and  B respectively;  the  segment  FE  is  a parallel  to  the  axis  of 
symmetry  and  located  at  a distance  of  120  from  the  latter.  As  the  main  equation  we  use  the  equation  of 
continuity 


d_ 

dx 


(pvvr) 


^-(pv,.r)=0 

or 


(3.1) 


where  vx  and  vr  are  the  velocity  components  in  the  x and  r directions  respectively.  The  flow 
velocity  and  liquid  density  are  scaled  by  the  free  stream  velocity  V0  and  density  p0  respectively.  In  the 
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compressible  case  equation  (3.1)  is  essentially  nonlinear,  and  we  will  apply  the  finite-difference  method 
to  solve  it.  The  density- velocity  relationship  is  taken  from  the  Bernoulli  equation 


P = 


) ,n  = 7A5 


(3.2) 


In  order  to  make  the  fonnulation  of  the  boundary  conditions  at  the  cavity  surface  more 
convenient,  it  makes  sense  to  map  the  computational  domain  in  x,  r coordinates  onto  a unit  square  in  rj 
coordinates,  conformity  of  the  points  is  shown  in  Fig.  5.  After  mapping  equation  (3.1)  takes  the  form 
[18] 
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( prUc  ] 

a 

( prV'  ) 
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dx 
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dx  dr  dr  dx 


(3.3) 


where  If  and  V are  the  contravariant  velocity  components  in  the  g and  l]  directions,  and  J is  the 
Jacobian  of  the  mapping.  We  consider  the  flow  to  be  potential  and  write  the  velocity  components  in  the 
following  fonn  [18] 


Uc 


= A, 


3<I> 


+ A 


3<1> 


3 dn  ’ 


vc 


= A 


ao 

¥ 


+ A-, 


ao 

a?7 


rd£,  ' 

2 

'3£ " 

2 

? ^2  — 

f a/?  Y (dr]  Y 

ydx  j 

+ 

Kdr  J 

1 ^ 
F 

1 ^ 

3 dx  dx  dr  dr 


where  <J>  is  the  flow  velocity  potential  and  Au  A2,  A 3 are  metric  coefficients.  After  mapping  onto 
£/?  plane  equation  (3.2)  takes  the  fonn 


P = 


(n  - 1 )M 


, H dll 


(3.4) 


The  boundary  condition  is  the  following:  the  tangential  velocity  is  constant  on  the  cavity 
boundary  MN,  the  impermeability  condition  is  imposed  on  the  segments  CA,  AM,  NB  and  BD,  the 
potential  on  the  external  boundary  is  equal  to  the  free  stream  potential. 
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Thus,  in  order  to  determine  the  cavitation  flow  it  is  necessary  to  solve  equation  (3.3)  onto  a unit 
square  with  the  given  boundary  conditions.  The  dependence  of  the  density  on  the  velocity  (3.4)  is  added 
to  equation  (3.3).  The  entire  computational  domain  (that  is  the  unit  square  in  the  rj  plane)  is  divided  by 
a network  with  steps  AE,  and  Arj  along  c,,  and  rj  directions. 

The  non-linear  system  of  equations  obtained  as  a result  of  the  discretization  was  solved  using  the 
iteration  method  with  approximate  factorization  [18].  Before  proceeding  to  the  calculations  the  first 
approximation  for  the  cavity  shape  between  the  points  M and  N was  preassigned  and  the  computational 
network  corresponding  with  the  given  cavity  was  constructed  using  an  algebraic  mapping  technique  [18]. 
The  velocity  on  the  cavity  surface  Vk  was  determined  on  each  iteration  stage  using  the  values  of  the 
potentials  @N,  <t>M  at  the  points  N and  M from  the  correlation 


Vk={<S>N-<S>M)lsN 


,v_v  is  the  length  of  the  arc  measured  from  the  disk  edge  M to  the  point  N on  the  closing  disk.  The 
potential  at  the  point  on  the  cavity  surface  corresponding  to  s is  given  by  the  expression 


®{s)=Vks  + ®M 


(3.5) 


where  5 is  the  arc  length  measured  from  the  point  M to  the  given  point  on  the  cavity  surface. 
This  expression  was  used  for  specifying  the  boundary  conditions  on  the  segment  MN. 

The  iteration  procedure  turned  out  to  be  convergent;  the  discrepancy  diminished,  and  the 
potential  increment  vanished  as  the  number  of  the  iteration  cycles  increased.  However,  in  general,  the 
obtained  solution  did  not  satisfy  the  impermeability  condition  on  the  cavity  surface.  In  order  to  satisfy 
the  impermeability  condition  the  cavity  shape  was  varied  using  the  following  way  [1 1],  The  normal  flow 
velocity  on  the  cavity  surface  was  represented  in  the  form  - dO/dn = d(p/dn + dx/dn  where  <p  is  the 
perturbed  velocity  potential,  x is  the  coordinate  of  a point  on  the  cavity  surface,  n is  the  normal  to  the 
surface,  and  dx/dn=dr/ds.  The  cavity  shape  must  satisfy  the  impenneability  condition  dO/dn  =0  or 


dr  _ d(p 
ds  dn 


(3.6) 


In  [11]  it  is  shown  that  for  subsonic  flow  in  view  of  cavity  symmetry  about  the  mid-section  the 
iteration  process  of  solving  the  system  of  two  first  order  differential  equations  converges  precisely  to 
relation  (3.6).  The  system  has  the  form 


dr 

ds 


d(p  dx 
dn  ’ ds 


1- 


' dr  ' 

yds  , 


(3.7) 


According  to  equations  (3.7)  the  front  half  of  the  cavity  had  to  be  corrected. 

The  numerical  analysis  was  performed  for  the  constant  cavity  length  Lk= 199.96  that  corresponds 
with  the  cavitation  number  0.02  for  incompressible  fluid.  The  results  of  the  numerical  calculations  were 
checked  on  the  satisfaction  of  the  mass  and  momentum  conservation  laws  and  were  compared  with  the 
results  obtained  by  the  slender  body  theory.  The  numerical  calculations  show  that  in  subsonic  flow  the 
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cavity  shape  is  close  to  an  ellipsoid  of  revolution  and  the  dependence  a(X,  M)  is  described  by  the  the 
second  expression  (1.5)  with  sufficient  accuracy.  Moreover,  the  test  calculations  were  performed  to 
compare  with  results  [15]  for  incompressible  fluid  (M=  0,  0=0.02)  and  results  [10]  (M=  0.5  0=0.16).  The 
comparison  has  shown  a good  agreement.  For  subsonic  flow  the  cavity  shape  in  the  vicinity  of  the  disk  is 
shown  in  Fig.  4 in  lecture  #1.  The  continuous  curve  represents  the  calculations  [15]  for  incompressible 
fluid.  The  author’s  numerical  results  are  shown  by  points  for  different  Mach  numbers.  In  Fig.  5 in 
lecture  #1  we  compare  the  cavity  profiles  in  compressible  and  incompressible  fluid  for  the  same 
cavitation  number  0-0.0235.  The  continuous  curve  corresponds  to  the  cavity  profile  in  compressible 
fluid  calculated  by  the  numerical  method  for  M=  0.8;  the  broken  curve  corresponds  to  the  cavity  in 
incompressible  fluid  calculated  by  Logvinovich  formula  [14]. 


4.  The  application  of  the  numerical  methods  for  calculation  of  supercavities  in 
supersonic  flow 

Let  us  consider  axisymmetric  supersonic  cavitation  flow  past  a disk  in  accordance  with  the 
Riabouchinsky  scheme.  The  proposed  flow  scheme  is  presented  in  Fig.  5.  In  case  of  supersonic  flow  the 
shock  waves  appear  (see  Fig.  5,  curve  1 is  the  forward  shock  formed  upstream  of  the  disk,  and  curve  2 is 
the  trailing  shock  emanating  from  the  closing  disk  edge).  The  forward  shock  is  the  most  intense  and 
nearly  normal  in  the  vicinity  of  the  symmetry  axis. 

However,  the  shocks  appearance  does  not  break  the  condition  that  the  flow  is  potential.  In  water 
over  a wide  range  of  pressures  a shock  adiabat  agrees  with  the  static  one  expressed  by  the  Tait  equation 
[19],  The  special  investigation  on  the  analysis  of  the  shock  equations  was  earned  out  [20],  The  analysis 
has  shown  that  we  can  assume  that  the  normal  shock  in  water  is  isentropic  and  the  flow  is  potential  if 
Mach  number  does  not  exceed  the  value  2.2.  In  [20]  it  is  shown  that  we  can  determine  the  density  from 
the  Bernoulli  equation  (3.2)  in  the  whole  area  of  the  supersonic  flow  (at  the  known  flow  velocity). 

Thus,  the  main  equations  for  the  calculation  of  the  supersonic  cavitation  flow  are  equations  (3.3) 
and  (3.4)  as  in  case  of  subsonic  flow.  We  will  apply  the  finite-difference  method  to  solve  these  equations 
[12].  As  distinct  from  the  subsonic  flow,  in  the  supersonic  flow  the  boundary  conditions  on  the  external 
boundary  are  the  following:  the  potential  is  equal  to  the  free  stream  potential  only  on  the  part  of  the 
external  boundary  where  the  normal  velocity  directed  along  the  external  normal  is  less  than  the  velocity 
of  sound.  On  the  part  of  the  external  boundary  where  the  normal  velocity  is  more  than  the  velocity  of 
sound  the  potential  is  determined  by  interpolation  of  values  in  the  nearest  nodes  of  the  computational 
network.  As  distinct  from  the  subsonic  flow,  we  apply  the  artificial  viscosity  so  that  the  difference 
scheme  will  be  stable  in  the  supersonic  area  [18,21].  The  artificial  viscosity  is  introduced  as  a 
modification  of  density  expression  [12].  The  modification  is  equivalent  to  upstream  displacement  of  the 
computation  domain  nodes. 

The  non-linear  system  of  equations  obtained  as  a result  of  the  discretization  was  solved  using  the 
iteration  method  with  approximate  factorization  [12]  (as  in  case  of  subsonic  flow).  The  iteration 
procedure  turned  out  to  be  convergent;  the  discrepancy  diminished,  and  the  potential  increment  vanished 
as  the  number  of  the  iteration  cycles  increased.  However,  in  general,  the  obtained  solution  did  not  satisfy 
the  impermeability  condition  on  the  cavity  surface.  In  case  of  subsonic  flow  the  impermeability 
condition  was  satisfied  by  varying  the  cavity  shape  with  the  help  of  equations  (3.7)  which  are  first  order 
differential  equations,  and  since  at  subsonic  flow  velocities  the  cavity  is  symmetric  with  respect  to  the 
mid-section,  a single  boundary  condition  sufficed  and  only  the  front  half  of  the  cavity  had  to  be 
corrected.  However,  at  supersonic  flow  velocities  the  cavity  is  asymmetric  and  its  shape  must  be 
corrected  using  second  order  differential  equations  subject  to  the  boundary  conditions 


5 = 0:  r = 1,  x = 0 
s = sN\  r = 1,  x = Lk 
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The  cavity  shape  was  corrected  from  the  disk-cavitator  edge  towards  the  closing  disk  edge.  The 
system  of  differential  equations  was  derived  by  differentiating  (3.7)  with  respect  to  5 


d2r  _ d d(p  | d2 x _ drd'r  ^ dr  V 
ds  ds  , dn  ) ds 2 ds  ds2  , ds  , 


(4.1) 


Note  that  the  second  derivative  d2r/dx2  has  a singularity  near  the  disk  edge,  an  accurate 
determination  of  the  second  derivatives  in  (4.1)  presents  considerable  difficulties  in  this  area.  The 
calculation  of  the  second  derivatives  near  the  disk  edge  was  avoided  by  using  equations  (4.1)  for  2 <x<Lk- 
2,  and  in  the  vicinity  of  the  cavitator  and  closing  disk  the  cavity  shape  was  corrected  using  equations 
(3.7)  (as  in  subsonic  flow).  As  a result  of  the  numerical  solution  the  obtained  cavity  shape  satisfied  both 
the  constant  velocity  condition  (3.5)  and  the  impermeability  condition  (3.6). 

The  numerical  analysis  was  performed  for  the  constant  cavity  length  I4-=199.96  that  corresponds 
with  the  cavitation  number  0.02  for  incompressible  fluid.  Calculations  of  supercavities  downstream  of  a 
disk  within  the  range  of  Mach  numbers  \<M<\A  have  been  performed.  As  a result  of  the  analysis  for 
different  Mach  numbers  the  following  characteristics  were  defined  - the  cavitation  number,  the  cavitation 
drag  ratio,  the  mid-section  radius,  the  aspect  ratio,  the  cavity  shape,  and  the  shock’s  distance  from  a disk. 

The  results  of  the  numerical  calculations  of  the  supercavities  in  supersonic  flow  satisfy  the  mass 
and  momentum  conservation  laws.  For  the  supersonic  flow  the  cavity  shape  defined  by  the  numerical 
calculation  is  close  to  an  ellipsoid  of  revolution.  It  agrees  with  the  author’s  results  obtained  with  the  help 
of  the  slender  body  theory  (see  chapter  2).  However,  the  results  on  asymmetry  of  the  cavity  shape  (mid- 
section displacement  relative  to  the  middle  of  the  cavity  length)  do  not  agree.  While  the  numerical 
calculation  yields  a displacement  directed  downstream,  in  chapter  2 (the  first  expression  (2.15))  an 
opposite  displacement  was  obtained.  The  behavior  of  the  dependence  of  a on  A and  M is  also  different: 
for  constant  Lk  the  calculated  aspect  ratio  and  a increase  with  the  increase  in  the  Mach  number  (as  in 
subsonic  flow),  whereas  the  second  expression  (2.15)  yields  the  opposite.  The  discrepancy  is  due  to  the 
fact  that  the  search  for  the  second  approximation  to  the  cavity  shape  and  the  dependence  a (X,  M)  carried 
out  in  chapter  2 overstepped  the  limits  of  the  slender  body  theory  for  supersonic  flow.  In  fact,  the  first 
approximation  (an  ellipsoid  of  revolution)  is  not  a slender  body  in  the  front  part  where  the  shock  arises. 
The  shock  is  not  described  by  the  small  perturbation  theory.  This  is  the  second  peculiarity  of  the  slender 
body  theory  application  to  the  supersonic  cavitation  flows  that  should  be  noted. 

The  cavity  profile  in  the  vicinity  of  the  disk  is  shown  in  Fig.  6. 


Fig.  6 


The  front  of  the  cavity 
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The  continuous  curve  represents  the  calculation  for  incompressible  fluid  [15].  The  author’s 
numerical  results  for  supersonic  flow  are  shown  by  points  for  different  Mach  numbers.  It  can  be  seen 
that  in  supersonic  flow  the  cavity  in  initial  regions  is  narrower  than  in  incompressible  flow.  Narrowing 
of  the  cavity  is  related  to  the  Prandtl-Meyer  stream  at  the  disk  edge. 
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Fig.  7 Cavity  profile 

In  Fig.  7 we  compare  the  cavity  profile  in  supersonic  flow  and  incompressible  fluid  for  the  same 
cavitation  number  <7=0.0308.  The  continuous  curve  corresponds  to  the  numerical  method  for  M=  1.2,  the 
broken  curve  corresponds  to  the  cavity  in  incompressible  fluid.  The  maximum  asymmetry  to  the  cavity 
shape  was  observed  for  the  Mach  number  equal  to  1.1  (this  case  is  illustrated  in  Fig.  8) 
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Fig.  8 Cavity  profile 


In  Fig.  8 curve  1 corresponds  to  the  numerical  calculation,  curve  2 corresponds  to  an  ellipsoid  of 
revolution  (the  origin  is  on  the  disk  edge).  From  Fig.  8 it  is  evident  that  in  the  first  approximation  the 
cavity  profile  is  close  to  an  ellipsoid  of  revolution  in  supersonic  flow. 
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The  positions  of  the  forward  shock  and  the  sonic  line  (where  the  local  Mach  number  equals  1 ) 
were  determined  by  numerical  calculation.  It  was  found  that  in  case  of  supersonic  cavitation  flow  past  a 
disk  the  distance  between  the  forward  shock  and  the  disk  surface  is  much  greater  that  in  case  of 
continuous  supersonic  air  flow  [12]. 


PI  = 1.1 


Fig.  9 Position  of  the  forward  shock  for  Af=l.l.  1 is  the  forward  shock, 
2 is  the  sonic  line  and  3 is  the  cavity  front  profile. 


As  an  example,  in  Fig.  9 the  position  of  the  forward  shock  is  shown  for  the  Mach  number  equal 
to  1.1.  Curve  1 corresponds  to  the  forward  shock;  curve  2 corresponds  to  the  sonic  line:  curve  3 
corresponds  to  the  cavity  front  profile.  The  Mach  number  dependence  of  the  cavitation  drag  coefficient 
of  a disk  at  zero  cavitation  number  Cx„  is  shown  in  Fig.  10  [12]. 
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Fig.  10  Dependence  of  Cxo  on  M 


Based  on  the  numerical  calculations  we  can  conclude  that  the  passing  through  the  speed  of  sound 
does  not  result  in  a substantial  variation  in  the  cavity  shape  as  compared  with  subsonic  flow.  In 
agreement  with  the  slender  body  theory,  the  cavity  shape  in  supersonic  flow  is  close  to  an  ellipsoid  of 
revolution. 
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5.  Shocks  in  a supersonic  water  stream 


We  will  consider  an  oblique  shock  in  a supersonic  water  flow,  the  double  line  in  Fig.  11,  where  /3 
is  the  angle  between  the  direction  of  the  undisturbed  velocity  V0  and  the  shock,  6 is  the  angle  of 
deviation  of  the  flow  behind  the  shock,  vno  and  vw  are  the  normal  and  tangential  components  of  the 
velocity  ahead  of  the  shock,  and  vni  and  vTi  are  the  velocity  components  behind  the  shock. 


Fig.  11  Oblique  shock  and  diagram  of  supersonic  flow  past  a cone  OA 


The  basic  relations  for  the  shock,  namely,  the  mass  and  momentum  conservation  equations  can 
be  written  as  follows  [22] 


Pov»  = PlV„i 

(5.1) 

py„o +p0= piv„2, + p, 

(5.2) 

> 

II 

t-T 

> 

(5.3) 

where  P0,  p0  are  the  pressure  and  the  density  ahead  of  the  shock,  and  subscript  1 denotes  the 
quantities  behind  the  shock. 

The  equation  of  the  shock  adiabatic  curve  must  be  added  to  equations  (5.1)-(5.3).  It  is  well 
known  [19]  that  at  the  pressures  lower  than  3- 1 03  MPa  the  shock  and  the  static  adiabatic  curves  coincide 
and  can  be  expressed  by  the  Tait  equation 


P.-P o=B 


Pj_ 

yPo  j 


-1 


5=  — f^o  n = 7A5 


(5.4) 


where  tfo  is  the  speed  of  sound  in  the  free  stream  ahead  of  the  shock. 
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If  on  the  shock  front  the  pressure  exceeds  the  value  3- 103  MPa,  the  equation  of  the  shock 
adiabatic  curve  can  be  written  in  the  form  [19] 


Pi  - P0  = d 


/ _ \ 


\P o J 


, d = 4\6MPa,  k = 6.29  (5.5) 

Behind  the  shock  the  speed  of  sound  (a,2  = dP / dp)  can  be  determined  from  the  static  adiabat 


(5.4) 


a 


2 

1 


fp.  r 


(5.6) 


The  expressions  for  determining  the  flow  velocities  ahead  of  the  shock  Vo  and  behind  it  V\  have 
the  form 


+ v. 


TO  9 


+ v 


2 

rl 


The  dependences  of  the  ratio  of  densities  po/pi  on  the  free  stream  Mach  number  (M=  Vo/ao) 
come  from  equations  (5.1)-(5.3),  (5.4)  and  (5.5)  and  can  be  written  as  follows 


M2  sin2  (3 
M2  sin2  f3 
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( ~ V 

Pi 
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— 1 
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j _ 

vPo  y 

8 P'  l 
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f1--! 
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(5.7) 


We  apply  the  first  expression  (5.7)  when  on  the  shock  front  the  pressure  does  not  exceed  the 
value  3-103  MPa  (Msinp< 2.2),  the  second  expression  (5.7)  corresponds  to  the  pressure  that  exceeds  this 
value. 

Let  us  consider  a particular  case  that  corresponds  to  the  normal  shock  (j3=  90°).  For  this  case  the 
following  correlations  are  satisfied 


r0 


0,  Vn 


Vx  = vBl , sin  13  = 1 


The  quantity  of  perturbed  velocity  u behind  the  shock  is  determined  from  the  expression 


u = V0-V, 


(5.8) 
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The  following  is  the  way  how  we  can  calculate  the  flow  parameters  behind  the  normal  shock. 

The  value  of  the  pressure  on  the  shock  front  P\  is  given,  then  from  (5.4)  or  (5.5)  the  ratio  of  densities 

pi/ Po  is  calculated.  From  expressions  (5.7)  the  Mach  number  and  velocity  Vo  are  determined  , from  (5.1) 

the  value  of  velocity  behind  the  shock  V\  is  found.  Then  from  equations  (5.6)  and  (5.8)  the  speed  of 

sound  behind  the  shock  a\  and  perturbed  velocity  u are  determined.  With  the  help  of  formulas  (5.1), 

(5.2),  (5.4)-(5.8)  the  author  has  calculated  the  flow  parameters  behind  the  normal  shock  on  the  pressure 

P i interval  490MPa<Pi<7840  MPa  [20].  The  results  of  this  calculation  have  been  compared  with  the 

results  [23]  obtained  from  the  mass,  momentum  and  energy  conservation  equations  which  were  solved  by 

the  method  of  successive  approximations.  The  flow  parameters  calculated  by  the  author  (Vo,  u,  a\),  and 

/ / 

determined  in  [23]  at  the  same  values  P\  ( V0  ,u',  a,  ) are  presented  in  the  table. 


Table 


Pi, 

MPa 

Vo,  m/s 

Vo',  m/s 

u,  m/s 

u',  m/s 

a\,  m/s 

a \,  m/s 

490 

1967 

1975 

249 

251 

2221 

2230 

980 

2310 

2335 

424 

426 

2734 

2755 

2630 

568 

567 

3175 

2880 

694 

689 

3535 

2450 

3033 

3110 

808 

798 

3798 

3855 

2940 

3223 

3320 

912 

898 

4075 

4140 

3430 

3394 

3510 

1011 

990 

4343 

4405 

3920 

3550 

3690 

1104 

1075 

4605 

4650 

4900 

3836 

4020 

1277 

1235 

5088 

5100 

5880 

4095 

4325 

1436 

1380 

5527 

5505 

6860 

4333 

4610 

1583 

1510 

5931 

5880 

7840 

4554 

4885 

1722 

1625 

6309 

6240 

From  the  table  it  is  evident  that  the  results  obtained  by  the  author  is  close  to  the  results  obtained 
in  [23].  Consequently  we  can  use  expressions  (5.4)  and  (5.5)  for  calculation  of  the  shock  adiabatic  curve, 
ft  should  be  noted  that  T.  Nishiyama  in  [9]  applied  the  equations  identical  to  the  relations  for  a normal 
shock  in  air  for  calculation  of  a normal  shock  in  water.  In  [24]  it  is  shown  that  the  equations  in  [9]  are 
not  applicable  for  calculating  the  shocks  in  water  and  they  contradict  with  the  well-known  experimental 
and  theoretical  results. 


16-20 


Let  us  continue  the  consideration  of  the  oblique  shock.  Expressions  for  the  flow  turning  angle  9 
and  the  Mach  number  M,  = Vl  / at  behind  the  shock  can  be  determined  from  equations  (5.1)-(5.3)  and 
(5.6)  and  the  geometric  relations  in  Fig.  1 1. 


tgjfi-e) 

tgp 


M sin  p 
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sin 


(P-9) 


Pcl 

P. 


\(»+l)/2 


(5.9) 


The  flow  turning  angle  9 was  calculated  at  Mach  numbers  on  the  interval  1.2<M<5  [25].  The 
calculation  was  performed  as  follows.  For  the  chosen  value  of  M the  interval  of  shock  angle  p (from 
Mach  angle  equal  to  arcsin  1/M  to  n/2)  was  divided  into  200  parts;  for  the  chosen  values  of  M and  p the 
ratio  of  densities  pjp\  was  determined  from  (5.7)  by  half-division  method.  The  value  of  pjp\  was 
substituted  in  the  first  equation  (5.9)  and  the  value  of  9 was  determined.  The  Mach  number  M\  was 
determined  from  the  second  equation  (5.9).  In  Fig.  12  we  have  reproduced  the  calculation  results  in  the 
form  of  graphs  of  9 versus  p at  various  M (continuous  curves  1-8,  the  Mach  numbers  equal  1.2;  1.4;  1.6; 
2;  2.5;  3;  4;  5 respectively). 


Fig.12  Dependence  of  6 on  /3;  curves  1-8  correspond  to-M=1.2;  1.4;  1.6;  2;  2.5;  3;  4 and  5; 
curve  9 and  10  correspond  to  M)=l  and -0=9*  respectively 

The  chain  curve  9 separates  the  domains  of  supersonic  flow  behind  the  shock  (M\>  1)  and 
subsonic  flow,  the  domains  of  supersonic  flow  are  located  to  the  left  curve  9.  The  broken  curve  10 
corresponds  to  the  maximum  (critical)  angles  0*.  The  domain  to  the  right  of  curve  10  corresponds  to  a 
detached  shock.  If  the  apex  angle  of  the  wedge  exceeds  the  value  of  20*  then  the  llow  around  this  body 
has  the  curved  detached  shock.  In  Fig.  13  curve  1 represents  the  graph  of  the  critical  angle  0*  as  a 
function  of  M for  a wedge. 
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Fig.13  Dependence  of  0*  on  M;  curves  1 and  2 correspond  to  flow 
past  a wedge  and  cone  respectively 


6.  Conical  flows 

Let  us  consider  the  axisymmetric  conical  flow  formed  behind  a shock  (Fig.  1 1)  [25];  the  line  OA 
is  the  profile  of  the  cone  and  a is  its  apex  half-angle.  The  spherical  coordinate  system  (p,  r with  the 
origin  that  lies  at  the  apex  of  the  cone  and  with  x-axis  directed  along  the  symmetry  axis  is  introduced. 
The  angle  <p  is  measured  from  the  x-axis,  /--coordinate  equals  the  radius-vector  drawn  from  the  origin,  OB 
is  the  intermediate  conical  surface,  Vr  is  the  component  of  the  velocity  directed  along  the  radius-vector, 
F(p  is  the  component  of  the  velocity  directed  along  the  perpendicular  to  the  radius-vector. 

The  basic  condition  of  the  conical  flow  concludes  that  the  flow  parameters  must  be  constant 
along  the  straight  lines  drawn  from  the  apex  of  the  cone.  Consequently,  any  partial  derivative  of  the  flow 
parameters  with  respect  to  r - coordinate  equals  zero.  The  equations  that  define  the  kinematic  parameters 
of  the  conical  flow  have  been  obt  ained  from  the  continuity  and  Euler  equations  [22] 


dVv  _2a2Vr+a%ctg<p-VrV2  dV^ 
d(p  V2-a2  ’ d(p 


where  a is  the  speed  of  sound  at  the  given  point  of  the  flow. 

For  the  velocity  components  Vvi  and  Vrl  the  boundary  conditions  on  the  shock  are  written  as 

follows 


Vt  = V0  cos {3,  q>  = p 
Pi 


(6.2) 


The  boundary  conditions  on  the  cone  surface  have  the  form 


r,=o,  <p  = oc 


(6.3) 
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where  Vcon  is  the  quantity  of  the  velocity  on  the  cone  surface. 

The  speed  of  sound  in  the  region  behind  the  shock  is  determined  from  the  Bernoulli  equation 


2 


n - 1 


+ 


v?h- 


n - 1 


(6.4) 


The  dependence  of  the  density  on  the  speed  of  sound  has  the  form 


P = 


P^ 


f \2/(n— l) 

a 


(6.5) 


The  pressure  in  regions 
known  value  of  the  density 


behind  the  shock  is  calculated  from  the  static 


adiabatic  curve  at  the 


P = 


P\a\ 


-1 


(6.6) 


The  conical  flow  behind  the  shock  has  been  calculated  in  the  following  way.  The  system  of 
equations  (6.1),  the  boundary  conditions  (6.2)  and  (6.3)  and  equation  (6.4)  were  reduced  to  the 
dimensionless  form.  The  quantity  of  the  velocity  was  normalized  by  the  free  stream  speed  of  sound  a0 
and  we  obtained  the  Mach  number  M as  the  dimensionless  parameter.  The  Mach  number  was  specified 
on  the  interval  \<M< 5.  For  the  chosen  M the  interval  of  shock  angle  /3  (from  Mach  number  to  n/2)  was 
divided  into  200  parts.  Then  for  the  chosen  values  of  M and  /3  from  (5.7)  and  (6.  l)-(6.4)  the  apex  half- 
angle of  the  cone  a was  calculated.  The  system  of  differential  equations  (6.1)  was  solved  by  the 
numerical  Runge-Kutta  method.  As  the  boundary  condition  we  applied  condition  (6.2)  and  the  step  of 
angle  change  \(p  was  equal  to  0.25°.  The  condition  (6.3)  was  used  to  determine  the  end  of  counting.  The 
segment  [q>„  <p,iJ  was  determined  on  which  the  value  of  V<t  reversed  sign.  Then  on  this  segment  the 
value  of  a that  corresponds  to  the  condition  Vv=0  was  found.  As  a result  of  these  calculations  for  the 
fixed  values  of  M we  have  obtained  the  dependences  a(j5),  Mcon(P),  Vco„fl3)  (Mco„=Vco„/aco„,  where  acon  is 
the  speed  of  sound  on  the  cone  surface). 

From  the  dependences  afp)  the  maximum  angles  a corresponding  to  the  critical  cone  angles  9 * 
were  determined  for  all  values  of  M.  If  the  cone  apex  half-angle  a exceeds  the  value  of  9*  then  the 
supersonic  flow  around  this  body  has  the  curved  detached  shock  and  the  conical  flow  behind  the  shock  is 
absent.  In  Fig.  13  curve  2 represents  the  graph  of  9*  as  a function  of  M.  A comparison  of  the  results  in 
Fig.  13  with  the  results  for  air  [22]  shows  that  in  supersonic  water  flow  past  a wedge  or  cone  the  critical 
angles  are  significantly  less  than  those  in  an  air  flow  past  the  same  bodies. 

For  the  visual  representation  the  found  values  of  [i  a , Vco„  and  Mco„  were  treated  and  the 
dependences  Vco„(M)  and  Mcon(M)  were  obtained  for  the  fixed  values  of  a.  The  treatment  was 

performed  by  interpolation.  The  values  of  cone  angles  a were  chosen  equal  to  5°,  10°,  15°,  20°,  25°  and 
30°.  Using  the  value  of  the  velocity  on  the  cone  surface  Vcon  the  pressure  coefficient  Cp  was  calculated 
from  equations  (6.4),  (6.5)  and  (6.6). 

In  Fig.  14  we  have  plotted  graphs  of  the  angle  /3  of  inclination  of  the  shock  as  a function  of  M at 
cone  angles  5°<ct<300  (curves  1-6). 
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Fig.  14  Dependence  of  j8  on  M;  curves  1-6  correspond  to  cc=  5°,  10°,  15°,  20°,  25°  and  30°;  curve  7 
corresponds  to  detachment  of  the  shock  from  the  apex  of  the  cone;  curve  8 corresponds 
to  Mcon= 1;  curve  9 corresponds  to  the  Mach  line  ( sinp=l/M ) 

The  broken  curve  7 corresponds  to  detachment  of  the  bow  shock  from  the  apex  of  the  cone.  The 
chain  curve  8 separates  the  domains  of  supersonic  {Mcon>  1)  and  subsonic  flows  on  the  cone,  the 
supersonic  domain  is  located  below  curve  8.  The  broken  curve  9 corresponds  to  the  Mach  angle 
(, sinp=\/M ).  From  Fig.  14  it  is  evident  that  the  bow  shock  for  the  cone  o=5°  is  Mach  shock,  it  is  close  to 
the  Mach  line.  We  can  assume  that  the  flow  around  this  cone  is  described  by  the  small  perturbation 
theory. 


Fig.  15  Dependence  of  Cp  on  M;  curves  1-6  correspond  to  a =5°,  10°,  15°,  20°,  25°  and  30°;  curve 
7 corresponds  to  the  detachment  of  the  shock  from  the  apex  of  the  cone;  curve  8 

corresponds  to  Mcon= 1 
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In  Fig.  15  we  have  plotted  graphs  of  the  pressure  coefficient  Cp  as  a function  of  M for  various 
cone  angles  (curves  1-6).  The  broken  curve  7 and  the  chain  curve  8 have  the  same  meaning  as  those  in 
Fig.  14;  the  supersonic  domain  is  located  below  curve  8.  In  case  of  supersonic  flow  around  the  cone  the 
drag  coefficient  C,  equals  the  pressure  coefficient  Cp  at  oK)  (see  chapter  2). 


Fig.  16  Dependence  of  Cxo  on  M for  the  cone  a=  5°;  curve  1 corresponds  to  numerical 
calculation;  curve  2 corresponds  to  calculation  from  formula  (2.10) 


Let  us  compare  the  drag  coefficient  of  the  cones  at  <7=0  (Cxo=Cp)  obtained  in  this  chapter  with  the 
ones  calculated  based  on  the  slender  body  theory  (2.10).  In  Fig.  16  and  Fig.  17  the  dependences  CX0(M) 
are  shown  for  the  cones  a=  5°  and  a=  1 0°  respectively.  Curve  1 corresponds  to  the  values  of  Cx„  obtained 
in  this  chapter,  curve  2 corresponds  to  the  values  of  Cxo  calculated  by  formula  (2.10). 


Fig.  17  Dependence  of  Cxo  on  M for  the  cone  a=  10°;  curve  1 corresponds  to  numerical 
calculation;  curve  2 corresponds  to  calculation  from  formula  (2.10) 

From  Fig.  16  it  follows  that  for  the  cone  a=5°  the  results  of  the  numerical  calculation  obtained  by 
formulas  (6.1  )-(6.6)  are  close  to  the  results  obtained  by  formula  (2.10)  at  the  Mach  numbers  l.l<A/<3. 
This  conclusion  agrees  with  Fig.  14  which  shows  that  the  bow  shock  for  the  cone  0=5°  is  close  to  the 
Mach  line. 
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From  Fig.  17  we  can  see  the  significant  discrepancy  between  the  values  of  Cxo  calculated  by 
formulas  (6.1)-(6.6)  and  the  ones  calculated  by  formula  (2.10).  We  cannot  consider  the  cone  ce=10°  to  be 
a slender  body  in  supersonic  water  flow.  The  numerical  calculation  shows  (Fig.  14,  15)  that  at  the  Mach 
numbers  on  the  interval  \<M<\2  this  cone  has  a detached  shock  and  the  velocity  on  the  cone  surface 
exceeds  the  speed  of  sound  if  the  Mach  number  is  greater  than  1.46.  This  flow  is  not  described  by  the 
small  perturbation  theory.  The  third  peculiarity  of  the  slender  body  theory  application  to  the  supersonic 
cavitation  flow  is  restriction  of  its  results.  For  the  drag  coefficient  of  thin  cone-cavitators  this  theory 
provides  the  exact  values  if  the  cone  apex  half-angle  does  not  exceed  5°. 


7.  The  comparison  of  the  theoretical  results  with  the  experimental  data 

In  the  experimental  works  [26,  27,  28]  the  cavitation  flows  formed  at  the  high  speed  body  motion 
in  water  have  been  investigated  on  the  Mach  number  interval  0.1<M<0.93  (the  body  has  a disk-cavitator). 
From  the  treatment  of  the  experimental  data  we  can  conclude  that  in  subsonic  flow  on  the  whole  range  of 
the  Mach  numbers  the  cavity  is  symmetrical  with  respect  to  the  mid-section  and  the  cavity  profile  is  close 
to  an  ellipsoid  of  revolution.  This  conclusion  agrees  with  the  author’s  results  obtained  with  the  help  of 
the  slender  body  theory  and  the  numerical  methods.  In  Fig.  18  the  continuous  curve  represents  an 
ellipsoid  of  revolution;  the  points  correspond  to  the  experimental  data  [28]  obtained  on  the  Mach  number 
interval  0.3<M<0.93;  the  cavitation  numbers  have  the  order  10‘4h-10’. 


Fig.  18  Cavity  profile;  the  continuous  curve  corresponds  to  an  ellipsoid  of  revolution,  the  points 

correspond  to  the  experimental  data  [28] 


In  [26,  27]  the  experimental  values  of  the  cavity  dimensions  obtained  on  the  Mach  number 
interval  0. 1 <M<0.7  and  the  cavitation  number  interval  10"4<CT<10 2 have  been  compared  with  the 
asymptotic  dependences  at  (7— >0  [29].  The  asymptotic  formulas  for  the  cavity  mid-section  radius  and  the 
cavity  length  scaled  by  the  disk  radius  are  written  as  follows  [29] 


R 


k 


cx  lu- 
ff 


(7.1) 


In  Fig.  19  the  continuous  curve  represents  the  asymptotic  dependence  of  the  mid-section  radius 
on  the  cavitation  number  (7.1);  points  1 correspond  to  the  experimental  data  [26,  27]  (the  Mach  numbers 
are  marked  above  the  abscissa  axis);  points  2 correspond  to  the  author’s  calculations  [20]. 
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Fig.  19  Dependence  of  Rk  on  <T,  the  continuous  curve  corresponds  to  dependence  (7.1);  points  1 
correspond  to  the  experimental  data  [26,  27];  points  2 correspond  to  calculations  [20] 


In  Fig.  20  the  continuous  curve  represents  the  asymptotic  dependence  of  the  cavity  length  on  the 
cavitation  number  (7.1);  points  l correspond  to  the  experimental  data  [26,  27];  points  2 correspond  to  the 
author’s  calculations  [20], 
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Fig.  20  Dependence  of  Lk  on  <?,  the  continuous  curve  corresponds  to  dependence  (7.1);  points  1 
correspond  to  the  experimental  data  [26,  27];  points  2 correspond  to  calculations  [20] 


From  the  examination  of  Fig.  18-20  we  can  conclude  that  the  theoretical  results  agree  with  the 
experimental  data. 

In  the  United  States  and  in  the  Institute  of  Hydromechanics  (Ukraine)  the  tests  are  being  carried 
out  on  the  supersonic  body  motion  in  water.  The  quantitative  experimental  data  are  not  presented, 
however,  in  1998  a photograph  of  supersonic  cavitation  flow  with  a detached  shock  was  published  [30]. 
Also,  the  announcement  was  made  in  [30]  that  in  the  United  States  on  July  17,  1997  the  underwater 
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speed  record  was  set  (for  the  first  time  in  water  the  speed  of  a body  exceeds  the  speed  of  sound).  The 
contour  of  the  shock  on  the  photograph  shows  that  the  Mach  number  is  equal  to  1.1  (or  the  speed  of  a 
body  exceeds  1600m/s).  In  general,  the  image  of  the  flow  on  the  photograph  corresponds  to  the  author’s 
theoretical  result  for  the  Mach  number  equal  to  1 . 1 (Fig.  9). 


Conclusions 


The  high  speed  cavitation  flow  in  water  has  some  peculiarities.  The  cavity  has  remarkable 
characteristics;  it  promotes  not  only  the  motion  in  water  with  small  drag  but  also  screens  the  body  from 
shocks.  On  the  cavity  surface  the  constant  velocity  and  the  constant  pressure  conditions  are  satisfied, 
consequently  the  shocks  are  absent.  The  bow  shock  forms  before  the  cavitator  or  on  its  apex  (in  case  of 
thin  cones),  in  the  domain  of  cavity  closure  the  trailing  shock  forms.  For  this  scheme  of  flow  the  sharp 
change  of  the  flow  parameters  does  not  take  place  on  the  transonic  velocity  range  (0.7<M<1.2).  As 
distinct  from  water,  in  air  on  the  transonic  velocity  range  the  shocks  form  on  the  body  surface  and 
aerodynamic  characteristics  of  body  change  considerably.  In  water  on  the  transonic  velocity  range  the 
parameters  of  the  cavitation  flow  change  smoothly  as  the  Mach  number  increases.  The  passing  though 
the  speed  of  sound  does  not  result  in  a substantial  variation  in  the  cavity  shape  as  compared  with 
subsonic  flow,  in  spite  of  the  small  asymmetry  of  the  shape  about  the  mid-section  the  cavity  shape  in 
supersonic  flow  is  close  to  an  ellipsoid  of  revolution  as  it  was  before. 

Water  is  rather  law-compressible  fluid  as  compared  with  air.  In  water  many  physical  effects  tied 
up  with  compressibility  essentially  differ  from  those  in  air.  We  can  assume  that  the  normal  shock  in 
water  is  isentropic  and  the  flows  are  potential  at  the  Mach  numbers  on  the  interval  \<M<2.2.  At  the 
same  Mach  numbers  the  ratio  of  densities  p)/po  on  the  shock  front  in  water  is  far  less  than  that  in  air.  As 
a result  the  distance  between  the  forward  shock  and  the  cavitator  surface  in  water  is  much  greater  than  in 
case  of  continuous  supersonic  air  flow.  In  supersonic  water  flow  past  a wedge  or  a cone  the  critical 
angles  are  significantly  less  than  those  in  an  air  flow  past  the  same  bodies.  The  flow  turning  angle  for  the 
Prandtl-Meyer  stream  in  water  is  far  less  than  that  in  air  at  the  same  Mach  numbers.  All  the  factors 
mentioned  above  (the  small  losses  of  mechanical  energy  in  the  shocks,  the  considerable  distance  between 
the  forward  shock  and  the  cavitator  surface,  the  small  deflection  of  flow  in  the  stream  about  the  cavitator 
edge)  weakly  prevent  from  the  cavity  expansion.  In  supersonic  flow  at  the  Mach  numbers  on  the  interval 
1<M<1.4  we  do  not  observe  the  considerable  narrowing  and  change  of  the  cavity  shape  as  compared  with 
subsonic  flow.  The  cavity  expands  in  accordance  with  the  law  of  conservation  of  energy  in  a liquid. 

The  slender  body  theory  has  a limited  application  for  the  investigation  of  the  supersonic 
cavitation  flows.  First  of  all,  in  case  of  the  thin  cone-cavitators  the  solution  obtained  for  the  whole  area 
of  the  flow  is  not  valid  for  the  flow  past  the  cone  edge.  In  the  supersonic  flow  the  Prandtl-Meyer  stream 
forms  at  the  edge.  The  stream  must  turn  at  the  angle  that  is  dependent  upon  the  pressure  within  the  cavity 
or  the  cavitation  number.  The  incline  of  free  stream  line  is  less  than  the  cone  apex  half-angle,  ft  is 
impossible  to  satisfy  this  condition  within  the  limits  of  the  slender  body  theory.  Secondly,  for  the  outer 
solution  the  search  for  the  second  approximation  to  the  cavity  shape  went  beyond  the  limits  of  the  slender 
body  theory.  In  fact,  the  first  approximation  (an  ellipsoid  of  revolution)  is  not  a slender  body  in  the  front 
part  where  the  shock  forms.  The  shock  is  not  described  by  the  small  perturbation  theory.  Thirdly,  for  the 
drag  coefficient  of  thin  cone-cavitators  this  theory  provides  the  exact  results  if  the  cone  apex  half-angle 
does  not  exceed  5°. 

As  distinct  from  the  supersonic  flow  the  slender  body  theory  is  more  applicable  for  investigation 
on  subsonic  water  flow.  For  example,  the  author  has  determined  the  cavitation  parameters  in  subsonic 
flow  around  thin  cones  with  apex  half-angle  equal  to  5°,  10°  and  15°.  These  results  satisfy  the 
momentum  conservation  law  and  for  the  particular  case  of  incompressible  fluid  (M=  0)  agree  with  the 
numerical  calculations.  For  the  case  when  the  parameter  of  the  cavity  thinness  is  far  less  than  the 
parameter  of  the  cavitator  thinness  the  outer  solution  was  found  for  subsonic  flow.  The  comparison 
between  the  results  of  outer  solution  and  the  results  of  numerical  calculation  has  shown  good  agreement. 


16-28 


References 


1.  Gurevich  M.  I.,  “Half-Body  with  Finite  Drag  in  Subsonic  Flow”,  Tr.TsAGI,  #653,  1947,  pp  1-12. 
(in  Russian) 

2.  Yakimov  Yu.  L.,  “Asymptotic  Laws  of  Degeneration  of  Thin  Cavity  Shapes”,  Izv.  Akad.  Nauk 
SSSR,  Mekh.  Zhidk.  Gaza,  #3,  1981,  pp  3-10.  (reprinted  in  the  USA) 

3.  Yakimov  Yu.  L.,  “Thin  Cavitation  Cavity  in  a Compressible  Fluid”,  in  Problems  of  Contemporary 
Mechanics,  Pt.l,  Moscow  University  Press,  Moscow,  1983,  pp  63-73.  (in  Russian) 

4.  Vasin  A.  D.,  “Thin  Axisymmetric  Cavities  in  Subsonic  Compressible  Flow”,  Izv.  Akad.  Nauk 
SSSR,  Mekh.  Zhidk.  Gaza,  #5,  1987,  pp  174-177.  (reprinted  in  the  USA) 

5.  Vasin  A.  D.,  “Thin  Axisymmetric  Cavities  in  Supersonic  Flow”,  Izv.  Akad.  Nauk  SSSR,  Mekh. 
Zhidk.  Gaza,  #1,  1989,  pp  179-181.  (reprinted  in  the  USA) 

6.  Serebryakov  V.  V.,  “Asymptotic  Solutions  of  Problems  of  Axisymmetrical  Supercavitation  Flow 
in  Slender  Body  Approximation”,  in  Hydrodynamics  of  High  Speeds,  Chuvash.  University  Press, 
Cheboksari,  1990,  pp  99-11 1.  (in  Russian) 

7.  Serebryakov  V.  V.,  “Asymptotic  Solutions  of  Axisymmetrical  Problems  of  Subsonic  and 
Supersonic  Separated  Water  Flows  with  Zero  Cavitation  Numbers”,  Dokl.  Akad.  Nauk  Ukraine, 
#9,  1992,  pp  66-71  (in  Russian) 

8.  Al’ev  G.  A.,  “Separated  Transonic  Water  Flow  past  a Circular  Cone”,  Izv.  Akad.  Nauk  SSSR, 
Mekh.  Zhidk.  Gaza,  #2,  1983,  pp  152-154.  (reprinted  in  the  USA) 

9.  Nishiyama  T.  and  Khan  O.  F.,  "Compressibility  Effects  on  Cavitation  in  High  Speed  Liquid  Flow 
(Second  report  - transonic  and  supersonic  liquid  flows)”,  Bull.  JSME,  v.  24,  #190,  1981,  pp  655- 
661. 

10.  Zigangareeva  L.  M.  and  Kiselev  O.  M.,  “Calculation  of  the  Compressible  Subsonic  Cavitation 
Flow  past  a Circular  Cone”,  Prikl.  Mat.  Mekh.,  v.58,  #4,  1994,  pp. 93-1 07.  (reprinted  in  the  USA) 

11.  Vasin  A.  D.,  “Calculation  of  Axisymmetric  Cavities  Downstream  of  a Disk  in  Subsonic 
Compressible  Fluid  Flow”,  Izv.  Rus.  Akad.  Nauk,  Mekh.  Zhidk.  Gaza,  #2,  1996,  pp  94-103. 
(reprinted  in  the  USA) 

12.  Vasin  A.  D.,  “Calculation  of  Axisymmetric  Cavities  Downstream  of  a Disk  in  Supersonic  Flow”, 
Izv.  Rus.  Akad.  Nauk,  Mekh.  Zhidk.  Gaza,  #4,  1997,  pp  54-62.  (reprinted  in  the  USA) 

13.  Vasin  A.  D.,  “Application  of  the  Slender  Body  Theory  to  Investigate  Developed  Axisymmetric 
Cavitation  Flows  in  a Subsonic  Stream  of  Compressible  Fluid”,  Applied  hydromechanics,  Kiev, 
v.2  (74),  #3,  2000,  pp  17-25.  (in  Russian) 

14.  Logvinovich  G.  V.,  “Hydrodynamics  of  Flows  with  Free  Boundaries”,  Naukova  Dumka,  Kiev, 
1 969,  2 1 5 p.  (in  Russian) 

15.  Guzevskii  L.  G.,  “Plane  and  Axisymmetric  Problems  of  Hydrodynamics  with  Free  Surface”, 
Dissertation  for  a doctor’s  degree,  Institute  of  Thermal  Physics,  Siberian  Branch  of  the  USSR 
Academy  of  Sciences,  Novosibirsk,  1987,  300  p.  (in  Russian). 

16.  Liepmann  H.  W.  and  Roshko  A.,  “Elements  of  Gasdynamics”,  John  Wiley  & Sons,  Inc.,  New 
York,  1957. 

17.  Ashley  H.  and  Landahl  M.  “Aerodynamics  of  Wings  and  Bodies”,  Addison-Wesley  Publishing 
Com.,  Inc.,  Massachusetts,  1966. 


18.  Fletcher  C.  A.  J.,  “Computational  Techniques  for  Fluid  Dynamics”,  Springer  Verlag,  Berlin  et  al., 
1988. 


16-29 


19.  Zamyshlyaev  B.  V.  and  Yakovlev  Yu.  S.,  “Dynamic  Loads  in  Underwater  Explosions”, 
Sudostroenie,  Leningrad,  1967,  387  p.  (in  Russian) 

20.  Vasin  A.  D.,  “Problems  of  Hydrodynamics  and  Hydroelasticity  at  High  Speed  Body  Motion  in 
Water”,  Dissertation  for  a doctor’s  degree,  TsAGI,  Moscow,  1999,  282  p.  (in  Russian) 

21.  Holst  T.  L.  and  Ballhaus  W.  F.,  “Fast,  Conservative  Schemes  for  the  Full  Potential  Equation 
Applied  to  Transonic  Flows”,  AIAA  Journal,  v.  17,  #2,  1979,  pp  145-152. 

22.  Loitsyanskii  L.  G.,  “Mechanics  of  Liquids  and  Gases”,  Pergamon  Press,  Oxford,  1966. 

23.  Cole  R.  H.,  “Underwater  Explosions”,  University  Press,  Princeton,  1948. 

24.  Vasin  A.  D.,  “High-Speed  Body  Motion  in  Compressible  Fluid”,  in  Proceedings  of  the  Scientific 
Meeting  on  High-Speed  Hydrodynamics  and  Supercavitation,  Laboratoire  des  Ecoulements 
Geophysiques  et  Industriels,  Grenoble,  France,  2000. 

25.  Vasin  A.  D.,  “Shocks  and  Conical  Flows  in  a Supersonic  Water  Stream”,  Izv.  Rus.  Akad.  Nauk, 
Mekh.  Zhidk.  Gaza,  #5,  1998,  pp  196-199  (reprinted  in  the  USA) 

26.  Savchenko  Yu.  N.,  Semenenko  V.  N.  and  Serebryakov  V.  V.,  “Experimental  Check  of 
Asymptotic  Formulas  for  Axisymmetric  Cavities  at  a—> 0”,  in  Problems  of  High-Speed 
Hydrodynamics,  Chuvash.  University  Press,  Cheboksari,  1993,  pp  117-122.  (in  Russian) 

27.  Savchenko  Yu.  N.,  Semenenko  V.  N.  and  Serebryakov  V.  V.,  “Experimental  Investigation  of 
Developed  Cavitating  Flows  at  Subsonic  Flow  Velocities”,  Dokl.  Akad.  Nauk  Ukraine,  #2,  1993, 
pp  64-69.  (in  Russian) 

28.  Savchenko  Yu.  N.,  “Investigation  of  High-Speed  Supercavitating  Underwater  Motion  of  Bodies”, 
in  “High-Speed  Body  Motion  in  Water”.  AGARD-R-827,  1998,  Reference  20. 

29.  Garabedian  P.  R.,  “The  Calculation  of  Axially  Symmetric  Cavities  and  Jets”,  Pacific  J.  Math.,  #6, 
1956,  pp  61 1-689. 

30.  Kirchner  I.  N.,  “Supercavitating  Projectile  Experiments  at  Supersonic  Speeds”,  Abstract,  in 
“High-Speed  Body  Motion  in  Water”,  AGARD-R-827,  1998,  Reference  35. 


